Raman response of Stage- 1 graphite intercalation compounds revisited 
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(Dated: March 3, 2013) 

We present a detailed in-situ Raman analysis of stage-1 KCg, CaCe, and LiCe graphite interca- 
lation compounds (GIC) to unravel their intrinsic finger print. Four main components were found 
between 1200 cm~^ and 1700 cm~^, and each of them were assigned to a corresponding vibrational 
mode. From a detailed line shape analysis of the intrinsic Fano-lines of the G- and D-line response we 
precisely determine the position {uiph), line width (Pph) and asymmetry (g) from each component. 
The comparison to the theoretical calculated line width and position of each component allow us to 
extract the electron-phonon coupling constant of these compounds. A coupling constant \ph < 0.06 
was obtained. This highlights that Raman active modes alone are not sufficient to explain the 
superconductivity within the electron-phonon coupling mechanism in CaCe and KCs. 



I. INTRODUCTION 



Sp^ hybridized carbon allotropes have unique struc- 
tural properties in different dimensions like graphite 
(3D), graphene (2D), single walled carbon nanotubes 
(ID), and fuUerene molecules ("quasi OD"). They have 
been widely studied due to their interesting electronic 
properties ranging from metallic, semimetallic, (zero gap) 
semiconducting to wide gap semiconducting and insulat- 
ing One unique possibility to tailor their electronic 
properties is by intercalation of alkali and alkaline-earth 
ions. These intercalation compounds are particular ap- 
pealing for their application in batteries and because of 
their superconducting phases. Superconductivity, as a 
result of alkali-metal intercalation, was first studied by 
Henning [5] in Graphite Intercalation Compounds (GIC) 
and further studies [1, Q- However, until 1981, the crit- 
ical temperature (Tc) in staqe-1 XCg GIC (X=K, Rb, 
and Cs) reported was low not higher than 0.135 K 
for CsCs, and between 0.39-0.55 K for KCg. This is sur- 
prising since GIC are BCS superconductors based on elec- 
tron phonon coupling owing an exceptionally high elec- 
tron phonon coupling constant up to A = 0.45 in the case 
of KCs and high phonon frequency of the optical modes 
Q. For instance, using A=0.45 and a phonon frequency 
of 1337 cm-i a BCS T^ of -5 K would be possible in KCg 
which is much higher than the observed Tc up to 0.55 K, 
and this can be related to a screened Coulomb pseudopo- 
tential of fj* = 0.14, which is on the lower bound with 
respect to CaCe 

In 1991, the discovery of fullerene intercalation com- 
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pounds, so-called fullerides, added a new family of or- 
ganic superconductors of type AaCgo (A=alkali-metal) 
jlOl4l2| . Compared to classical superconductors, and in 
contrast to GIC the Tc observed in fullerides is high rang- 
ing from 18 K for KgCeo, 28 K for RbgCeo, up to 39 K 
for CsaCgo [TTl - [l3| . Contrary to GIC, where the highest 
intercalation level represents the superconducting phase, 
for fullerides the superconducting phase is a line phase 
at half filling. Other stable fullerides AiCeo, A4C60 and 
AgCgo are either normal metals, Mott-Hubbard insula- 
tors or charge transfer insulators [l^. Similar to GIC 
the superconducting coupling mechanism was described 
within the framework of BCS theory involving an elec- 
tron phonon coupling to the intra-molecular modes of Ceo 
3 . Experimentally, most important for the coupling are 
the two low energy intra-molecular modes with Hg sym- 
metry [1^ [1^ , although theoretically the high frequency 
phonons have been predicted to play a significant role 

For GIC, the observation of superconductivity of CaCg 
in 2005 with a high Tc of 11.5 K [l?| triggered further 
research in the field and led to alternative explanations 
of the superconducting electron phonon coupling. For in- 
stance, Kim et al. attribute superconductivity in CaCg 
to the high-energy C modes [2^]. Hinks et al. [21| report 
that the low-energy modes of the intercalant were respon- 
sible for superconductivity inferred from specific heat 
analysis, while first principle calculations predicts equal 
coupling to both groups of phonons [1, [23 ■ Therefore, 
the exact contribution of the different coupling phonons 
still remain elusive. 

Raman spectroscopy became then an important tool 
to determine the exact contribution of each phonon, and 
it opened a route for revealing the coupling mechanism in 
superconducting fullerides and GICs. Hence, it serves as 
a key tool to analyze the electron phonon coupling con- 
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stant (A) from a renormalization of the optical response 
of the intra-molecular Cgo modes and of the graphitic 
G-Hne response. Recent Raman studies on the G-line 
response of different stage- 1 GIG reported the assign- 
ment of the electron phonon coupling (EPC) induced line 
width 7^^*^ to the 1510 cm"! mode [IS-IiB], which has 
been explained b y th e inclusion of non adiabatic phonon 
calculations (23. |26|[. However, the intrinsic G-line re- 
sponse in heavily doped graphite compounds are still elu- 
sive because of the influence of defects and laser induced 
deintercalation, as recently reported using a micro Ra- 
man analysis for GaGg 24!] and for KCg single crystals 

In this contribution we report a detailed study of the 
D- and G-lines in KCs, GaGg, and LiCg GIG, in order to 
unravel their intrinsic phonon components and its rela- 
tion to the electron phonon coupling constant responsible 
for superconductivity. From the analysis of the optical 
phonons observed, we assign their role in the supercon- 
ductivity coupling mechanism in comparison with previ- 
ous results of electron doped GIG. 



II. EXPERIMENTAL, AND MEASUREMENT 
DETAILS 

The synthesis of KGs was performed in-situ under high 
vacuum (^^4x10^^ mbar) conditions in a quartz tube 
with natural graphite flake single crystals from differ- 
ent sources, and a potassium ingot with 99.95% purity 
(Aldrich) for the intercalation. Potassium was evapo- 
rated until golden crystals were obtained. This phase can 
be directly assigned to stage-1 KGg phase from a com- 
parison of the Raman response with previous combined 
Raman and XRD results p7l. [28|. GaGg, and LiGe were 
prepared in a sealed ampoule by using a procedure de- 
scribed elsewhere [1^ . Highly oriented pyrolytic graphite 
(HOPG) flakes were degassed and used for lithium and 
calcium intercalation for 10 days under He atmosphere 
(ca. 0.5 atm). The ampoule was then opened in the 
glove box and gold colored product was extracted from 
the melt. Powder x-ray diffraction measurements were 
carried out using a Stadi-P diffractometer (GuK^) to con- 
firm the intercalation stage in GaGe and LiGg. For the 
Raman analysis every GIG was kept in vacuum (^4x10^'* 
mbar) in order to avoid de-intercalation due to expo- 
sure to air. The Raman analysis, was performed with a 
HORIBA LabRam at room temperature, with a 568 nm 
wave length, and 0.25 mW of laser power. Every spec- 
trum were acquired under the same conditions in a range 
from 500 cm~^ up to 2500 cm"-'^ and the line positions 
were calibrated by gauge lamps. 

III. RESULTS AND DISCUSSION 

In the Raman response of staqe-1 GIG eight optical 
vibrational modes are present |28| in the following irre- 
ducible representation: 
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FIG. 1. a) Optical modes of graphite, b) Raman spectra 
from KCg, CaCg, and LiCe taken with 568 nm laser at room 
temperature and low laser power of 0.25 mW. 



The Eo , and the Eo vibrational modes are Raman 
active, and the A2u, and Eiu belong to infra-red active 
modes P, H^l . There are some other modes in graphite 
which are forbidden in perfect graphite and only become 
active in the presence of disorder like the mode with Aig 
symmetry. In Fig. [T]a), the optical modes of graphite are 
depicted. Previous Raman studies in GIG have confirmed 
the presence of the E2g mode around 1600 cm"-'^, the 
{c-axis mode) around 500 cm^^, and the absence of 
the Aig 0,111]. The c-axis mode has being attributed 
to an out-of-plane G motion in graphite f23]. This mode 
correspond to the M point of the graphene Brillouin zone, 
and it becomes Raman active when high intercalation 
levels are achieved. In agreement with the literature, we 
observe (as shown in in Fig. [l]b), that the c-axis mode 
is present solely in KGg around '^560 cm^^. Surprisingly 
and in agreement with previous studies neither in CaGg, 
nor in LiGg this mode is observed [H, [3l| . 

Regarding the G-line response all these previous stud- 
ies reported one G-line which has a strong Fano line shape 
due to the coupling and the interference with the conduc- 
tion electrons. Taking a closer look on the lineshape of 
the G-line response in Fig. [l]b), one can easily see that 
more than one component is present, and a detailed line 
shape analysis is needed in order to unravel their intrinsic 
response and related electron phonon coupling of these 
stage-1 GIG. The line-shape analysis of the G-line is dis- 
cussed in detail below. 
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FIG. 2. D- and G-line analysis for stage-1 GIG, and Raman response of their laser induced deintercalated phases. The four 
components which can be identified in the G-line shape are: Aig mode between 1250 and 1350 cm~^, £292 mode ~1510 cm~^, 
E291 mode at ~1547 cm"^ , and stage-2 G-mode ~1560 cm~^. In the upper panel a) we can observe that KGs exhibit a strong 
contribution from the £2^2 with a broad Fano behavior, which is the finger print for the intrinsic line of stage-1 compound [27| . 
In the lower panel b) we present the same crystals analyzed in the upper panel but de-intercalated. We can clearly observe the 
decrease of the £292 mode, concomitant to a strong increase of the G-mode assigned to the XC24 graphitic face ^^1600 cm~^. 



A. Analysis of the intrinsic G-line response of 
Stage-1 GIG 

The structure of the intercalation stages in graphitic 
compounds, has been studied and it is well understood 
from x-ray diffraction ;28j. However, the intrinsic Ra- 
man response of stage-1 GIC is still complicated by laser 
induced de-intercalation from a local heating of the sam- 
ple with different laser power densities [13]. In addition, 
other factors such as 3D intrinsic disorder of the crys- 
tal also strongly affect the Raman response in GIC. For 
example, a graphite single crystal doped to stage-1 will 
remain polycrystalline due to a non-homogeneous inter- 
calation. This will limit the achievable doping in these 
GIC [11,133. 

Hence, the previous experimental and theoretical re- 
sults on the Raman response of KCs and CaCe reported 
in the literature are not conclusive with respect to the 
G-line shape and position. In different studies a wide 
range of different G-line positions between ~1400 cm~^ 
and ~1600 cm^^ are reported: i.e. at ~1500 cm~^ Uj, 
between 1400 cmT^ and 1550 cm"! 1534 cm"! [gg, 
1547 cm-i [lOl, 1420 cmT^ and 1582 cm-i_pl|. In more 
recent experiments for calcium GIC |24l |25| , potassium 
doped graphene and graphite [13, l32l|. and later in Li- 
graphite [31| , the strongest G-line phonon response is ob- 
served around 1510 cm~^ when the sample has the best 
quality (lowest defect content) and highest intercalation. 

In Fig. [5] a) the D- to G-band region of pristine stage- 
1 intercalation compounds with K, Ca, and Li is de- 



picted and clearly shows the presence of shoulders in the 
response, which indicate different components. Never- 
theless, in order to compare with the previous studies 
[26I [30I [33} we first conducted a line-shape analysis of the 
G-line by using a single Breit Wigner Fano (BWF) func- 
tion. This yields parameters which are in good agreement 
to those results, and confirms that our samples have the 
same high quality of a true stage-1 compound. This is 
further supported by the fact that the G-line assigned to 
stage-2 compounds around 1600 cm~^ is only increased 
upon e.g. laser induced de-intercalation (see Fig. [2]b). 

In a second step a detailed and accurate analysis of 
the line-shape in the D- to G-band region of these GIC 
was conducted using four components. The assignment 
of each component to the Aig(D), -£232' ^2gi modes, and 
the G-line of the stage-2 compound is explained in the 
following. 

Regarding the line-shape, all components have been 
fitted using BWF functions of the form: 



I{w) - /o 



(1 



qr/2 ) 



l + (3f^)2 



A 



where ujph is the phonon frequency, F the line width or 
damping, q the asymmetry parameter and A an offset. 
For the first and fourth peak (D, and G), the asymme- 
try was q=10^ approaching a Lorentzian function, while 
the second and third (splitted G-line) have a pronounced 
Fano interference. In the analysis, in order to get com- 
parable results for each GIC, the same values of F, and 
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q were used to fit each respective component. The pa- 
rameters are summarized in Table I together with the 
calculated values from the adiabatic and non-adiabatic 
phonons from Ref. [2^ . 

The first mode observed in Fig. [2] a) between 1260 
and 1360 cm~^ has been previously attributed to parti- 
cle size effects and/or the presence of disorder [s^, HH]. It 
has been assigned to the Aig vibration, which is forbid- 
den in perfect graphite. Therefore, this mode is called 
D-line (intrinsic "defect mediated"), and it involves the 
contribution from the phonons near the K zone boundary 
with a Lorentzian line-shape. 

The second and third modes observed are assigned to 
the E2g graphitic mode of heavily doped graphene lay- 
ers (Fig. [2] a). Both components have a pronounced 
asymmetry and they are well described by a BWF line- 
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shape. We label the two modes as E2gi 
E2gi mode is located between 1528 cm^^ and 1585 cm 
and it is attributed to not homogeneous or incomplete 
intercalation in stage- 1 compounds [27[. The £^292 mode 
locates at 1510 cm~^ for KCg and CaCe, and 1546 cm~^ 
for LiCg. It has a clear and strong Fano behavior which is 
characteristic to the finger print of stage- 1 graphite inter- 
calation compounds ^24j, |27j| . When de- intercalation was 
induced in the samples, a decrease of these i?2g modes 
was remarkably observed (see Fig. [Hb). 

The fourth mode related to the G-line of their respec- 
tive stage-2 compound is observed at 1612 cm^^ for KCs, 
1600 cm~^ for LiCg and at 1560 cm~^ for CaCg. The sur- 
prising low frequency in the case of CaCg was also found 
in Ref. [2^ and explained as a de-intercalated phase in 
CaCg. As mentioned above, the increase of this fourth 
component is highlighted in the partly de-intercalated 
stage- 1 compounds in Fig. [2]b), and points towards a 
phase separation upon de-intercalation. 

B. Analysis of the Electron-Phonon Coupling 

The previous results are very important for the correct 
determination of the stage, and electron-phonon coupling 
constant Xph responsible for superconductivity within the 
BCS theory [M, |26„ l36] . This constant is directly related 
to the intrinsic G-line phonon frequency, and to the adia- 
batic (uja) and non-adiabatic {ujj^a) phonon frequencies. 
Saitta et al. [26.] have analyzed the EPC in many dif- 
ferent stage- 1 GIG from a difference in the experimental 
phonon frequency to the calculated phonon frequency in 
the adiabatic and non-adiabatic limit. In order to de- 
termine the electron phonon scattering renormalized line 
width j^P^ [H, [2^ we used: 



„EPC 



(1) 



We obtain 7^^*^ values for KGg, GaGg, and LiGg 
which are in very good agreement to our experimental 
Vph value obtained from our BWF fit. Table II. In Fig. 
ini we show the location of our 7^^*^ with respect to the 
expected linear tendency to Tph as predicted by Saitta 



TABLE I. Fit parameters to the four components of the D- 
and G-line in the Raman spectra of KCg, CaCe, and LiCe. 
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et al. [26j. It is important to notice that some com- 
ponents of the G-line in KGg, GaGg, and LiGg bring a 
^EPC _ -ftrhich means that they do not show the non- 
adiabatic effects for layered metals and therefore they do 
not contribute to the electron-phonon coupling constant 
Xph. In comparison to the experimental F^^^* and 7-^^"-' 
from Ref. 0, IH [11] (Fig. E]*), our resuhs using the 



E2c,2 mode are in better agreement to the linear trend 
expected for F k, 7^^*^. This confirms the importance 
of every optical mode in the range between the adiabatic 
and non-adiabatic frequency range (w^-wat^), and con- 
firms that the E2g2 component is the intrinsic stage- 1 
vibrational mode with the strongest non-adiabatic effect 
on the EPG. 

We now turn to a detailed analysis of the EPG constant 
\ph. Different values have been already reported and 
used to calculate the critical temperature of KGs, GaGg, 
and LiGg with values around 5 K, 11.5 K, and 0.9 K, 
respectively, in agreement with some experimental and 
theoretical studies [1, |40| . In order to extract Xph from 
the phonon line-width (F) and position {ijjph) from our 
Raman data we used [4l|: 



A F2 



2 M ujr,K vjr 



(2) 



where the electron-phonon coupling strength is given 
by Dexp ■■ 



A n2 



(3) 



and Auc is defined as the area of the graphene unit 
cell, M is the carbon atom mass, vp is the Fermi velocity, 
AF(3 is the Landau damping phonon decay rate given by 
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FIG. 3. Calculated 7^^*^ (Eq. [1} for different GIC as function 

•k) correspond to experimen- 
. The dashed line represents 
the approximation otTph ~ "f^^^ . The red dots show our cal- 
culated EPC, which are in better agreement to the expected 
approximation to the Tph values. 

ATq = Tph - T Graphite, and Fpjf has dimensionality of 
a force taking in consideration the lattice displacement 
along the corresponding optical phonon mode. By using 
Eq. m and the definition of = 4:{DI)f, and F^. = 
2{D'j^)p from Ref. [sl, |4l| we calculate the values for 
Xr,K for each phonon in the T-K branch observed in the 
G-line region as summarized in the right column of Table 
II. (Dp j^)p were taken from the DFTgga calculations 
in Graphite [s^ as they are closer to our electron-phonon 
coupling strength {Dexp)- 

By using the averaged electron-phonon coupling con- 
stant Xph = Ar + Xk, and the position ujph from the 
strongest optical mode in KCg, CaCe, and LiCe one 
can estimate the critical temperature Tc using McMil- 
lan's formula [42]. Taking our bjph values converted in 
to phonon temperature 8, /i* sa 0.14 from [4^, and Xph 
from the Raman analysis, we obtain Xph < 0.06 values, 
which are too low to explain superconductivity within 
EPC mechanism using these high-frequency Raman ac- 
tive modes. 

However, this is not a general behavior in interca- 
lation compounds. Electron-phonon studies in alkali- 
intercalated fuUerenes showed the possibility to attribute 
the strongest Xph contribution for superconductivity to 
the Hg(l) mode in A3C60 fullerides [H, [l^. More over, 
in agreement to the analysis reported by Yao et al. in 
Ref. [l6| our D^xp presented the same trend as the one 
observed in fullerides intercalation compounds. There- 
fore, we can confirm that the larger the value of 1/q, the 
weaker the coupling strength Dg^p in GIC and fullerides. 

On the other hand, in comparison to the EPC constant 



TABLE II. Electron-phonon coupling parameters from the 
G-line Raman analysis. The values of uiph, Tph, 'y'^'''^ are 
in cm~^ and they were extracted from the BWF analysis of 
the Raman spectrum. Dexp is the electron-phonon coupling 
strength from Eq. [3] in (eV/A). 
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phonon decay in dressed electron-hole pairs "y^^^ Ref. [2^ 
Optical branch assignment based in [si. |39|| . 
Electron-phonon coupling constant from Eq. [2] 



Xarpes reported using an analysis of the self energy re- 
sults in ARPES [1, HJ] , our Xph values are about a factor 
of 10-15 lower. Since, in the case of CaCg superconduc- 
tivity was confirmed at Tc=11.5 K, only the Xarpes 
would be sufficient to explain this high superconducting 
transition temperature. Hence, the low Xph proves that 
optical modes from the G-line in stage- 1 GIC are not 
sufhcient to explain Tc in the electron-phonon driven su- 
perconducting coupling mechanism and additional not 
optically active modes might play an important role. 

IV. CONCLUSIONS 

We have performed a detailed in-situ Raman study of 
the most common GIC (KCg, CaCg, and LiCg). We iden- 
tify four main peaks in the D- to G-band region, and all 
these Raman responses match the spread of different line 
shapes reported in the literatures so far. From an evalua- 
tion of the fine structure in the G-line response we assign 
each peak to their corresponding vibrational mode and 
phonon branch. 

We found the strongest Fano behavior of the G-line at 
1510 cm~^ in KCs and CaCg, not like in LiCg, which 
highlights the importance of this mode to the supercon- 
ductivity coupling mechanism within the BSC theory, 
and confirms the importance of this i?2g2 mode to non- 
adiabatic effects. By using this mode, we obtain a very 
good agreement to the theoretical predicted line-width 
ryEPC ^ especially for CaCe. 

Finally, we find a very small EPC Xph < 0.06 which is 
much too low to explain the high Tc in this graphite 
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intercalated compounds. This points out that, other 
phonons including acoustic modes and other electronic 
states might play an important role in explaining the su- 
perconducting pairing in GIC. 
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